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INTRODDCTION 


There  are  three  versions  of  the  finite  element  method:  the  h-version, 
the  reversion  and  the  h-p  version.  The  h-version  is  the  standard  one,  where 
the  degree  p  of  the  elements  is  fixed,  usually  on  low  level,  typically  p  = 

1,2,3  and  the  accuracy  is  achieved  by  properly  refining  the  mesh.  The  in¬ 

version,  in  contrast,  fixes  the  mesh  and  achieves  the  accuracy  by  increasing 
the  degrees  p  of  the  elements  uniformly  or  selectively.  The  h-p  version  is 
the  combination  of  both. 

The  standard  h-version  has  been  thoroughly  investigated,  and  many 
programs  are  available,  both  commercial  and  research  codes.  The  p-version  and 
h-p  versions  are  new  developments.  There  is  only  one  commercial  code,  the 
system  PROBE  (Noetic  Tech,  St.  Louis).  The  theoretical  aspects  have  been 
studied  only  recently.  The  first  theoretical  paper  appeared  in  1981  (see 
[7]).  See  also  [1],  [7],  [8],  [9J,  (10J,  [11].  For  the  numerical,  computa¬ 
tional  and  irapleraentational  aspects  we  refer  to  [2], 

In  [4]  it  has  been  shown  that  the  rate  of  convergence  is  an  optimal  one 
up  to  an  arbitrarily  small  e  >  0,  namely 

(1.1)  lien  ,  <  c(e)p'(k-i  )+Gnun  . 

H 1  HK 

In  the  case  when  the  solution  has  singular  behaviour  of  the  type 
u  ■  r01  ’  a  >  0,  ((r,9)  being  polar  coordinates)  and  the  vertex  of  the 

elements  is  at  the  origin,  then 

(1.2)  Hell  <  C(e)p~2a+e. 

H 

The  p-version  has  the  rate  of  convergence  which  is  twice  the  rate  of  the  h- 
version  with  uniform  mesh. 


1 


proof . 

In  [4],  only  the  case  when  essential  boundary  conditions  are  homogeneous 
was  addressed.  In  this  paper  we  deal  with  the  general  case.  The  ideas  and 
techniques  of  this  paper  differ  significantly  from  what  was  used  in  [4].  Sec¬ 
tion  2  addresses  the  preliminaries  and  basic  notions.  Section  3  deals  with 
approximation  properties  of  polynomials  on  a  square.  Section  4  analyzes  the 
rate  of  convergence  under  the  assumption  that  the  solution  does  not  have  sin¬ 
gular  behaviour  and  proves  (1.3)  for  homogeneous  and  nonhoraogeneous  essential 
boundary  conditions.  Section  5  deals  with  the  case  when  the  solution  has  a 
singular  behaviour  and  proves  (1.4).  Section  6  summarizes  the  results  and 
addresses  briefly  various  generalizations. 


v  ^  y-y  •£££ * v:svv:-sy  d 


2.  PRELIMINARIES 


2.1.  Notation 

2  2 
By  R  we  denote  the  usual  Euclidean  space  with  x  =  (x^,X2)  €  R  . 

2 

By  fl  c  R  we  denote  a  bounded  polygonal  domain  with  the  vertices  A^ ,  i  * 

M 

0,...,M,  A q  m  Aw,  and  the  boundary  r  a  £  T  where  r.  are  open  straight 

i-1 

lines  with  the  end  points  Ai-i*Af  The  internal  angle  of  and  is 

denoted  by  u>^,  i  m  0  <  <  2it.  We  mention  that  we  also  allow 

(i)^  *  2n  and  so  we  include  into  our  consideration  the  slit  domains  (uk  *  2 tt ) 

when  the  boundary  is  two  sided  (in  an  obvious  sense). 

Let  rD  -  l  r.  and  rN  -  T  -  rD  *  l  f.,  V  n  W  =  0.  We  will  call 
jeP  J  jew  J 

n  J  xj  j 

the  Dirichlet  boundary  and  the  Neumann  boundary.  Obviously, 

rD  u  rN  -  r. 

By  L2(ft)  *  H9(fl)  and  ,  k  >  0,  integer,  we  denote  the  standard 

Sobolev  spaces  (with  index  2).  Also,  Hp(fl)  m  Hk(fl)  fl  H^(fl)  where  H^(fl)  = 

{u  €  H^(fl)|u  “0  on  rD} .  For  k  >  0  not  an  integer,  we  define  Hk(fl), 

H0(fl)  as  the  usual  interpolation  space  (by  the  K-method,  see  [6]): 


h£+9(a)  -  (H£(fl),H£+1(n)) 

0,q 

If 

with  a*  2,  0  <  0  <  1,  £  +  0  ■  k.  For  k  >  1  we  define  (fl)  = 

H*C+^(Q)  fl  H^(fl).  In  [3]  we  have  shown  that  if  <  2tt  thus 

H^9(«)  -  (H^),H^+1(«))0jq. 


Later  we  will  also  use  q  *  <*>  and  will  explicitly  mention  this  case. 

We  will  also  deal  with  the  Sobolev  spaces  H^r^),  H^(I),  I  =  (a.b't 
which  are  defined  for  k  integer  in  the  analogous  way. 


The  spaces  Hk(ft)  ,  Hp(ft), 
inner  products  will  be  denoted  by 
For  <  >  0  we  let 


Hk(ri)  ,  etc., 
^  ’  )Hk(n)> 


are  Hilbert  spaces  and  their 
etc. 


(2.1)  R(<)  -  {(x1,x2)  |  jxj  <  < ,  |  x2 1  <  <} 

k.  k 

and  by  H  (R(k))  c  h  (R(k))  we  denote  the  space  of  all  periodic  functions 
PER 

with  period  2tc. 

By  Pp(fl)  ,  respectively  Tp(R(tc)),  we  denote  the  space  of  all  algebraic, 
respectively  trigonometric  (with  period  2x) ,  polynomials  of  degree  at 
most  p  in  each  variable  on  12,  respectively  R(k).  Analogously  we 
define  P  (rt),  P  (I),  (I  -  (a,b)). 


2.2.  The  model  problem  and 

its  properties 

We  will  consider  the 

following  model 

problem 

(2.2) 

-Au  +  u  *  F 

on  ft 

(2.3a) 

u  =  g  on 

rD 

(2.3b) 

3u  , 

- —  “  b  on 

3n 

rN. 

The  model  problem  (2.2)  (2.3)  is  a  classical  case  of  the  elliptic  equation 
problem  on  a  nonsraooth  domain.  The  structure  of  this  problem  is  well 
studied.  We  refer  here  to  [13]  and  the  survey  paper  [15]  where  the  relevant 
information  and  references  could  be  found. 

We  shall  assume  that  the  solution  of  (2.2),  (2.3)  C3n  be  written  in  the 


following  form: 


where 


ux  €  Hj(O) 


q  >  1 


u 2  €  Hk(£2)  ,  u  -  g  on  r°, 


k  > 


(2.5) 


[i] 


I  c 


[i] 


£=1 


YCil  a[i] 

log  rj  1  rt£  $[il(e.)XEil(ri)  €  H*(£2) 


with  aEil  >  0,  YEi5  >  0,  *Eil(0i)  and  x^t^)  are  C° 


(or  sufficiently  smooth)  functions,  XEiE(r^)  ■  1  for  0  <  ^  <  ■£, 

XE^(r^)  =*  0  for  r^  >  2pE^.  By  (r^,9^)  we  have  denoted  the  polar  coordi¬ 
nates  with  the  origin  at  the  vertex  of  the  polygon  12.  The  partition 


(2.5)  is  typical  for  the  regularity  of  the  solution  of  the  problem  (2.2) 

ti] 


(2.3).  The  functions  u^  describe  the  singular  behaviour  of  the  solution 
caused  by  the  corners  of  £2  or  by  the  abrupt  changes  of  boundary  conditions. 
Function  U2  relates  to  the  nonhoraogeneous  Dirichlet  conditions  on  rD  and 
uj  relates  to  the  solution  of  the  problem  with  the  homogeneous  Dirichlet 
conditions.  For  the  details  and  proofs  of  the  partition  (2.5)  we  refer  to 
[13]  [14]. 

So  far  we  assumed  that  £2  is  a  polygon  and  we  considered  only  the  model 
problem  (2.2)  (2.3).  In  Section  6  we  will  make  comments  about  more  general 


cases . 


2. 3.  The  p-verslon  of  the  finite  element  method 
N 

Let  £2  =  U  £2 .  where  £2^  are  (open)  triangles  or  parallelograms. 

i“  1  1 

(In  Section  6  we  will  comment  on  curvilinear  triangles  and  quadrilaterals.) 


We  shall  assume  that  £2 .  fl  £2 .  =*  0  for  i  f  j,  £2  H  £2 .  is  either  the  empty 


r 


y 


iiw»wiwwjiii.¥.wiw?wmT  ^’.v  w»  ■  j  ■>' 


set  or  an  entire  side  or  a  vertex  of  ftj.  and  ftj .  We  will  assume  that  all 


vertices  of  ft  are  the  vertices  of  some  ft^ .  ft^  will  be  called  elements. 


Denote 


S  -  {u  €  H1 (ft)  |  u  €  P  (ft.),  i  = 
p  pi 


sD  *  s  n  Hi(n). 

p  p  D 


Let  c  be  a  side  of  the  element  ft^  and  let  A.  ,A.  €  rD  be  the  end 

1  J  Jl  ^ 2 


point  of  7.*  We  define  g  ^  €  P  (7-): 

J  P  P  3 


(2.6a) 


Jj] 


g„  J (A.  )  =  g( A.  ), 


1,2 


(2.6b) 


/  (g^Vh'ds  =  /  g'h'ds 


7, 


for  all  h  €  P  (7-)  with  h(A^  )  *  0.  We  define  then 
P  J  Ji 


(2.6c) 


gp  •  {sp31k 


t3  f 


The  p-version  of  the  finite  element  method  consists  now  of  finding  up  €  S  , 


Up  =  gp  on  T  such  that 


(2.7) 


(u  ,v)  .  =  /  Fv  dft  +  /  bv  ds 


H  (ft)  ft 


hold  for  all  v  £  S  . 

P 


(2.7)  is  the  usual  definition  of  the  finite  element  solution  when 


replacing  g  by  gp  on  Tu  so  that  gp  is  the  trace  of  a  function  in  Sr 


Our  construction  of  gp  is  slightly  restrictive  because  we  assumed 


that  g  £  H 1  (7 j  )  and  not  g  €  Kk(7.j),  -j  <  k  <  1.  This  restriction  is  not 
important  in  practice.  It  could  be  avoided  at  the  expense  of  simplicity  of 
construction  of  gp. 

Remark.  If  ft^  is  a  parallelogram  then  P  (ft)  is  meant  as  the  set  of  poly¬ 
nomials  in  variables  which  are  parallel  to  the  sides  of  ft  . 


3.  APPROXIMATION  PROPERTIES  OF  Sp 

Let  i 

0  -  (-1,1)  x  (-1,1);  y., 

i  =  1,2, 3, 4  be  the  i 

sides 

of  Q  and 

y  5  be  the 

diagonal 

xl  =  x2  of  Q. 

LEMMA  3.1. 

Let  u  € 

Hk(Q).  Then  there  exists  a  sequence 

ZP  € 

Pp(Q),  p  = 

0,1,2, ... 

such  that 

(3.1)  for 

k  >  0, 

q  =  0, 1 ,  q  <  k: 

II  u-z  II  < 

P  Hq(0) 

Cp‘<k~‘>).iu>  „ 

Hk(Q) 

(3.2)  for 

k  >  j: 

ii  u-z  n  n  < 

p  h°(Yi) 

Cp  Hull 

H  (Q) 

i  =  1 , . .  . 

(3.3)  for 

k  > 

II  u-z  U  .  < 

P  H1^) 

r  — ( k— 3/ 2) 

Cp  BuP 

HK(Q) 

i  —  1 , . .  . 

(3.4)  for 

k  >  y 

and  any  x  €  Q: 

| (u-z  ) (x) |  < 

P 

C  p~  ( k- 1  )  II  n 

Hk(Q) 

The  constants  C  in 

(3.1)  -  (3.4)  depend  in  general  on  k 

but 

are  inde- 

pendent  of 

u  and 

P- 

Proof.  Let 

r0  >  1. 

Then  Q  c  R(tq) 

(see  (2.1)).  Since 

Q 

is  a  Lipschit 

domain,  there  exists 

an  extension  operator  T  mapping  H 

(0) 

into 

Hk(R(2rQ)) 

such  that 

(3.5a) 

Tu  =  0  on 

R(2rn)  -  R(~  rn) 

(3.5b)  ITu*  <  CIluB 

Hk(R(2r0))  Hk(Q) 

where  C  is  independent  of  u.  For  a  concrete  construction  of  T  we  refer 
to  [3], 

Let  $  be  the  one-to-one  mapping  of  R(-|-)  onto  R(2tq): 

R(2tq)  3  x  =  (xltx2)  =  $(5)  =  (2tq  sin  £  ,  2r^  sin  E,^) 

with  ,C2)  =  5  €  R(^)  . 

Further,  we  let 

R  =  <t_1[R(|  rQ)  ]  c  R(i) 

where  <t>  ^  denotes  the  inverse  mapping  of  $ . 

Let  v  =  Tu  and 

(3.6)  V(0  =  v(*(0). 

Because  of  (3.5a)  we  easily  see  that 

(3.7)  supp  V(5)  c 
In  addition  it  can  be  readily  seen  that 

(3.8a)  V(£)  is  a  periodic  function  with  period  2tt 

(3.8b)  » V  (C ) «  k  <  C II  v  II  .  <  C  n  u  II 

n  (R(ir ) )  Hk(R(y  rQ))  H  (Q) 

and  hence  V  f  HkER(R(ir)). 

(3.8c)  V(£)  is  a  symmetric  function  with  respect  to  the  lines  £  =  ±  — , 


-.'.si'-'.-. 


Let  us  expand  the  function  V  in  terms  of  its  Fourier  series 


(3.9) 


VU1'52>  ■  l  £  ajle 

J  =*-,»  £=*—00 


i(jC1+iC2) 


and  deno  te 


(3.10) 


V£1'C2)  '  ^  £  aJie 

y  j— p  Z  — p 


i(jC1+iC2) 


Obviously,  Vp  €  Tj(R(ir))  . 


We  have 


l/2  2k 


(3.11) 


Hk(R(*)) 


I  U  l2((i+i2^2)  ) 

j,£ 


where  ■  has  the  usual  meaning  of  equivalency.  (3.11)  yields  immediately 
for  0  <  q  <  k 


(3.12) 


IIV-V  II  <  Cp 

P  Hq(R(7T)) 


•(k_q)m 


H  (R(if)) 


(using  (3.8.b)) 


<  CD-(k_q)l 


Hk(Q) 


with  C  independent  of  u. 


i  *  1,...,4  be  the  sides  of  R(tt)  and  let  =  £->  he  one 


the  sides.  Then 


(3.13)  V(£L  ,l2)  -  Vp(C1  ,l2)  - 


|jl>p  ut<p  |jl<p  u!>p  Iil 


>p  |Z|>p 


i( jC1+£l2) 


where  for  |j|  >  p: 


(3.14) 


'k 


( 1  ]  ,  2 


<p 


a. 


2 


(.  ?  b*1 

U!<p 


(by  Schwarz  inequality) 

<  (  I  laiJ2(1+j2+£2)kK  l  a+j2**2)-*) 

i4<p  J£  i4<p 


A,( 


j  (i+j2)k 


■) 


CA.p 


-( 2k- 1 ) 


where  we  denote 


(3.15) 


Aj 


l  la  |  2d-Hj 2+£2)k. 


For  |j|  <  p: 


(3.16) 


|bi 


[2] 


‘  (uf 


>p 


Aj  u! 


>p 


(i+j2+*2) 


-k 


•  2  2  ~k 
<  CA.  /  ( 1+j  +x  )  dx 

J  P+1 


<  CAjI(k,p) 


where 


provided  that  k  >  V2  • 
Analogously 


(3.18) 


|  3  J |  <  CA.I(k,p)  <  CA  p-(2k_1) 

J  J  J 


provided  that  k  >  V2  •  Hence  for  i  *  1,2, 3, 4, 


(3.19) 


IIV‘VPH  0  . 

hu(y.) 


<  c[ 


,<‘'l  - 


lb'2ll  ♦ 


IjT<p  2  IjI>p  2 


b>31!  1 


<  Cp'(2k_l>  l  A 

j=-—  J 


Cp-(2k-U  I  I  |a,tU^2«2)k 

j  =—0O  £=—00  •* 


Cp  (2k_1)llVII2 

H  (R(*)) 


<  Cp  (2k_1)auii2 

Hk(Q) 


provided  that  k  >  V2  . 

Now  we  estimate  H  V— V  I!  ,  .We  have 

P  Hl(?t) 


(3.20)  3J-CV-V  >  -  l[  y  jbJ'U151 

dh  p  1jT>p  3 


■  jb121^'1  *  y  Jb!3’eiJ?1] 

<P  J  i j l>p  J 


and  analogously  as  in  (3.14),  for  |jj  >  p, 


12 


(3.21) 


|jb]13 12  < 

2  (1+j  r  J 


Consider  now  the  function 


(l+x‘) 


2xk  * 


f'(x)  =  - pr—r  (l-Oc-l)x2) 

(l+x2)k  1 


and  hence  for  k  >  3/2  and  x2  >  2  we  obviously  get  f'(x)  <  0.  Hence  f(x) 
is  a  decreasing  function  for  x  >  2  and  k  >  3/2. 

Hence  for  |j|  >  p,  p  >  1: 


(3.22) 


Mb"1!  <  A.  <  CA  p 

2  (l+(p+l)z)k  2  2 


-( 2k-3) 


provided  that  k  >  3/2.  For  j  <  p  we  get  analogously  as  in  (3.16), 


(3.23) 


I  21  | 2  <  CA  /  - ^2—  <  CA  p 

2  2  p+1  (l+jZ+xV  2 


-( 2k-3 ) 


provided  k  >  3/2. 

Finally,  for  j  >  p, 


2  od  2 

[3],2  „  n  »  t  _Td X 


(3.24)  |jbJJ|  <  CA  J  -  2  2~  k 

3  2  p+1  (1+j  +x  )k 


CA.[  /  ( - 2  \  - 2  X  2  ~k~  '  ^ dx  ^ 

2  p+1  (l+j2+x2r  (1+j  +xZ)K 


<  CA.  p 
J 


-( 2k- 3) 
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i? 


provided  that  k  >  3/2.  Hence  from  (3.20)  and  (3.19) 


i 


I 


(3.25) 


SV-V  fl 2 ,  <  Cp"(2k_3)Sull2, 

P  H1 (y±)  Hk(Q) 


Let  us  estimate  now  j  (V-V^)(?^  ,5^)  i  •  Because  V-Vp  €  HpER(R(ir)), 
i  without  loss 
we  get  for  k  >  3/2 


assume  without  loss  of  generality  that  (S^,^)  ^  Y^»  Using  (3.13)  and 


(3.26)  (  l  | bj 


1,1 


[1] 


<  t,  i  j2ib‘ur)(  t  -L) 

|j|>p  2  I j i >p  j 


>p 


,  _  -( 2k-3)  2  1  .  r  -2(k-l )  2 

<  Cp  Hull  ,  —  <  Cp  Hull 

Hfc(Q)  P  Hk(Q) 


Using  (3.13),  (3.16)  we  get  for  k  >  3/2 


(3.27) 


(  I  bU1)‘ 
|j|<P  J 


‘  (bl 


I .  [ 2 ] |i  .  r  -2(k-l )  2 

|b  |  Jp  <  Cp  null 

<p  2  H K(Q) 


Finally,  using  (3.13),  (3.24)  we  get  for  k  >  3/2 


(3.28) 


|bl.3]|)  <  cp'2^1^^ 

|j|>P  J  Hk(q) 


(.  I  r 


Combining  (3.26),  (3.27),  (3.28)  and  (3.13)  we  get  for  k  >  3/2 


(3.29) 


|(V-V  )(Cl,42)l  <  Cp-(k  1 )!ul  . 

P  Hk(Q ) 


we  can 
.22) 


Let  us  prove  now  (3.19),  (3.25)  for  Y5*  We  have 


(3.30) 


§ 

il 

3 


s 


(V 

§ 

t 

i 

a 


3 


l| 

h 

H* 


V(C,5)  -  Vp(5,5) 


hl>p  uk  bk  uk  bk  uk 


£  (C^  +  +  C^3heiq5 

L  q  q  q 


i(j+*K 


where 


(3.31a) 


(3.31b) 


(3.31c) 


Now  by  the  Schwarz  inequality 


j+*-q 

I  j|>ptUI<p 


j+£-q 

I  j  l<P,  UI>P 


j  +  £=*q 

| j  t>p, |£)>p 


(3.32) 


<  (  I 

j  +  £-p 

I  J|>P,UI<P 


2  7  9  k  9  9 

i»ul  {i+j 2+t2)  )(  I  0+jV)  ) 

J  j+£=q 

I  j  I >p  •  Ul  <p 


where 


(3.33) 


<  Aqllll(k,p,q) 


2  2  2  k 
-  I  Uul  (i+j2+i2) 

j+£-q 


Illl(k,p,q)  <  - N(p,q) 

(1+p  ) 


and  N(p,q)  is  the  number  of  terms  in  the  second  term  on  the  right  hand  side  of 
(3.32).  Obviously,  N  <  2p.  Hence 


(3.34) 


|c[,l|  <  CA  p-<2k-U 

q  q 


Analogously , 


(3.35) 


|c121|  <  CA  p-<2k-U 

q  q 


Finally , 


(3.36) 


.t3]  ,  2  k  (  \ 


|C  Jj!  <  A  (  I  il+W)  ) 

q  q  j+£-q 

I  j  l>p.  Ul>p 


<  CA  7  <  Cp"(2k  °A  . 

q  ii  2k  r  q 

^  p+1  x 


provided  that  k  >  V2  .  From  (3.33)  we  see  that 


(3.37) 


T  a  <  ciivir 

u  n 


<  C II  u  It 


H  (  R  (  7T  )  ) 


Hk(Q) 


(3.37)  together  with  (3.34),  (3.  35),  (3.  36)  and  (3.30)  yields 


(3.38) 


-(k-  V? ) 

nv-v  1  .  <  cP  '2  nun  , 

P  H°Cy^ )  Hk(Q) 


Sow  we  estimate  II V  —  V  II  ,  .  Using  (3.30)  we  have 

P  H2(y5) 


(3.39) 


77  (V-v  )(e  ,5) 

d  K  P 


l  iq(C  ^  1  ^  +  C^3')eiq^ 

L  q  q  q 


sssssssasaBoa 


2 


I  i(^Cq  1 1  I  <  Aq  l  - 2 — TT 

q  q  j  +  £=q  (1+j  +  £ ) 

I  j  1>p,  M  <p 


If  >  q^  then 


l  — hr  <  p  I 

j  +  £=q  (1+j  +£  )  j+£ 


!  j|>p,UI<p 


4  2  2  k. 

j  +  £=q  ( 1+j  ^+1) 


I  j l >p . I  £  I <P 


<  c  — 1 —TT  <  CP 

(1+p  ) 


-( 2k-3 ) 


If  p^  <  q^  then 


P2  <  q 2 


( j+£ ) 2  <  2(j2+£2) 


and  we  get 


y  - a -  < 

j+£=q  (l  +  jZ+lZr 

I  j  l>p.UI<p 


I  q  2  k  *  (2p)^  2k-2 ^ 

J+l7q,  f i+4-)  q 

I  j  l>p»  U!  <p  v  2 


<  Cp 


—  ( 2k— 3) 


provided  k  >  1. 


Hence , 


(3.40) 


|,Cll]|  <  CA  „-(2k-3). 

1  q  1  q 


Similarly , 


(3.41) 


hc[2i|  <  CA  p-<2k-3). 

'  q  1  q 


Finally , 


(3.42) 


I  <  ca  l  — Hi 

q  q  j+£=q  (1+j ~+£^ ) 

I  j  l>P»  UI>P 


ww: 


For  q4  <  p4  we  get 


(3.43) 


l 


j+£=q  (l+j2+£2)k 

I  j  l>p.  M>p 


<  p 


|j]>P  a+jV 


<  Cp 


-(2k- 3) 


If  q^  >  p^  then 


(3.44) 


l 


2  2  k 

j+£=*q  (1+j  +£  )K 

|  j  |>p,  |£|>p 


2  k 


<  Cp 


-( 2k-3) 


(i+a_j 


and  hence 


(3.45) 


I  qc 


[3] 


<  CA  p 

q 


-(2k- 3) 


Combining  (3.40),  (3.41),  (3.45)  and  (3.39),  (3.38),  we  get 


(3.46) 


IV-V  H  ,  i  Cp~(k~  3/2\uH  , 

P  h‘(?5)  Hk(Q) 


Because  of  (3.8c),  V^($  ^(x))  €  P^(Q).  Further,  $  is  a  regular  mapping 

of  R(tq)  on  Q,  (tq  <  y)  and  4>  (y  ^ )  p  Yj..  Hence  for  k  integer  (3.1) 
follows  immediately  from  (3.12),  (3.2)  from  (3.19)  and  (3.38),  (3.3)  from 
(3.25),  (3.46),  and  (3.4)  from  (3.29). 

There  is  a  known  one  dimensional  version  of  Lemma  (3.1).  For  the  proof 
see  e.g.  [9] 

LEMMA  3.2.  Let  I  *  (-1,1),  u  €  H^(I),  k  >  1.  Then  there  exists  a  polynomial 
Zp  €  P  p( I)  such  that 

(3.47)  u(±l)=Zp(±l) 


■ 

* 


4 


and  for  t  «  0,1: 


where  C  depends  on  k  but  not  on  A  and  p. 
For  the  proof  see  [5]  or  [9]. 
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4.  THE  CONVERGENCE  RATE  OF  THE  p- VERS ION:  THE  CASE  WHEN  u  €  H  (8) 

4.1.  The  case  k  >  3/2. 

In  this  section  we  will  analyze  the  rate  of  convergence  of  the  p-version 
in  the  case  that  the  exact  solution  u  €  Hk(8) ,  k  >  3/2. 

THEOREM  4.1.  Let  u  €  Hk(8),  k  >  3/2  be  the  solution  of  (2.2)  (2.3).  Then 
there  exists  Up  €  Sp+j,  Up  =  gp  on  P^  (see  2.6a,b,c)  such  that 

(4.1)  Hu-u  B  .  <  Cp“(k_1 } il  uH 

P  H  (8)  H  (8) 


where  C  depends  on  the  partition  of  8  and  on  k,  but  is  independent  of  u 
and  p . 

Proof .  Let  8^,  i  ■  1,...,N  be  the  elements  of  the  partition  of  8.  First 

let  us  construct  the  functions  z  ^  as  in  Lemma  3.1.  The  lemma  is  applicable 

P 

because  a  linear  transformation  maps  the  parallelogram  or  triangular  element  on 
a  square  or  on  a  right  angled  triangle,  preserving  the  polynomials.  Hence 


(4.2) 


i!  u-z  ^  a 


^a  <  Cp-(k_1)auii 

p  Hl(8.)  Hk(8.) 

i  i 


Using  (3.4)  we  can  assume  that  u  =  z  at  the  vertices  of  8.-  bv  adding  a 

p  1 

linear  (triangle)  or  bilinear  (parallelogram)  function  to  z^^. 

Let  now  y  =  ^  ^  and  A^,  A2  be  the  end  points  of  y .  Now  Zp  ^ 

¥  z^  on  y.  Denote  w..=z^-z^  on  y.  Then  because  z^  (A.)  = 

P  J*  P  P  P  1 

r  0 1 

u(A^)  =  zL  J(A^),  we  have  wj£^A^)  =  0.  Using  (3.2)  and  (3.3)  we  get  for  t 


(  l  a'* 


II  u 


-(k-  Vo -t) 


If  ftj  is  a  parallelogram,  then  we  can  assume  without  any  loss  of  generality 
that  flj  =  Q  =  (-1,1)  *  (-1,1)  and  y  -  {x  |  | x1 |  <  1 ,  x2  -  -1}. 

Let 


(4.4) 


* P(x2) 


•P(x2+1)  -2p 

2  -e 

l-e‘2p  ’ 


x9  i  I  =  (-1,1 


Then 


(4.5a) 


M  o  <  Cp_1/2 

P  H° (I  ) 


(4.5b) 


K  4>  f)  <  Cp^2 

P  H 1  (I  ) 


and 


(4.6) 


>p(-D  =  1,  “  o. 


Using  Lemma  3.3  with  k  =*  1  and  t  =  0,1,  there  is  a  ^(x-,)  f  Pp(I) 
such  that 


(4.7) 


— 2+t  r  4-2p(x2+l)  2  V2 

n ^  — d>  n  <  cp  [/  p  e  (i-x2)dx2] 


<  Cpt(/  ye  2pydy)  <  Cp1-  *. 
0 


Hence,  for 


C  •  „  =  w  <p  (x„)  €  P  (n. ) 

J £  ji  P  2  P  J 


we  get 


(4.8)  K.,«  ,  <  C  Iw  I  .  I*  I  0  +  lw  l  0  H  i  L  . 

j£  H1^)  J£  h\y)  P  H°(I)  j£  HU(y)  P  H 1  (I ) 


<  cf  Iw.  ,  ( B  <t>  R  _  +  1$  -i|>  «  n  ) 

1  j£  H1 (y)  P  H° (I )  P  P  H°(I) 


+  lw..(  .  (hi,  +  «4>  -i>  a  ,  )] 

]Ih°(y>  p  h‘ci)  p  p»'(D 


by  (4.3) 


<  c[p-(k-  V2\p-^„-1)  *  p-<k-l/2)Cpi/2+1)llul 


Hk(ft) 


<  Cp~^k  ^  B  ill 


Hk(ft) 


Repeating  this  process  for  all  four  sides  of  ftj  we  can  adjust  z ^  so  that 

the  continuity  across  y  Is  obtained  and  (4.2)  still  holds  on  ftj . 

So  far  we  assumed  that  ftj  was  a  rectangle,  how  let  ftj  be  a 
triangle.  Then  without  any  loss  of  generality  we  can  assume  that 


ft.  *  (x  |  0  <  Xj  <  1,  0  <  x,  <  Xj} 


Y  *  (x  |  0  <  Xj  <  1,  x2=*0) 


and  we  assume  that  Wj£(0)  *  Wj^(l)  = 
Le  t  now 


^U(xl’x2^  3  ^p(2x2-1  ^  (xl-x2)wj£(xl  )  +  ('1-xl  ^wj£(xj-x2^ 


Obviously,  is  a  polynomial  of  degree  at  most  p  +  1  in  each  variable 

which  vanishes  on  3ft.  -  y  (because  w.  (0)  =  w.  .(1)  “  0)  and  ^..(x.,0) 

J  J  *•  J  *  ^  J  ^  1 

*  w.  (x  )  .  Now  by  quite  similar  arguments  as  before  we  see  that 
J  *  ^  * 


'  >'.-v 


H  y 


h  (o) 


Adjusting  by  £  •  we  achieve  the  continuity  across  y.  The  continuity 

across  each  side  y  of  £2^  belonging  to  1^  can  be  similarly  obtained.  Thi 
completes  the  proof  for  g  =*  0  and  for  g  a  polynomial  of  degree  pq  <  p  on 
every  3£y  c 

If  g  is  general  we  can  proceed  quite  analogously.  Let  y  c  be  the 

side  of  the  element  Oj  with  the  end  points  ,  i  =  1,2.  By  Lemma  3.1  we 
have  for  t  =  0,1  and  k  >  3/2: 


[j]-uJ  <  CP_(k"t"1/2)™ul, 

P  HC(y)  n  (0) 


On  the  other  hand,  by  the  imbedding  theorem  we  have  for  s  >  1 


HS(y) 


lull  <  C II  u  II  i, 

H  (y)  "  ■ 


H  'z(0) 


Applying  Lemma  3.2  we  have  for  t  =  0,1: 


Hence , 


Jg-g^1!'  <  Cp_(k  <  Cp  (k_t)llull  y 

p  hc(y)  hk(y)  hk  /2(o) 


(4.9) 


jj_gu  j« 

p  p  HC(Y) 


.  r  -(k-t-  V2  ) ..  , 
<  Cp  **  Hu  l 


hV.) 


Zljl(A,)  -  g[J>(A.) 

pi  d  1 


so  that  we  can  construct  the  adjustment  of  zp  exactly  as  before, 
This  completes  the  proof. 


4.2.  The  case  1  <  k  <  3/2 


We  assumed  in  Section  4.1  that  u  €  Hk,  k  >  3/2.  Let  us  analyze  now  the 
general  case. 

THEOREM  4.2.  Let  u  €  Hk(12) ,  k  >  1,  be  the  solution  of  (2.2)  (2.3)  such  that 

u  =  u^  +  U2 
ki  k9 

uj  €  HjjUft),  u 2  €  H  *(0),  k2  >  3/2 

and  such  that  if  k^  <  3/2  then  12  is  a  Lipschitz  domain  (i.e.  <  2i). 

Then  there  exists  up  €  Sp+p  up  *  gp  on  (see  2.6a,b,c)  such  that 

(4.10)  llu-u  II  ,  <  Cp“^k-1^llull  ,  ,  k  =  min(k,,k9) 

^  H  1  (ft )  HkU) 

Proof .  Because  of  Theorem  4.1  we  can  assume  that  U2  =  0,  i.e.  g  =  0  on  rD 
and  that  1  <  k  <  2.  Let  us  assume  first  that  for  any  0  <  t  <  °»  we  can  write 


(4.11) 


v,  (t)  +  v9(t) 


where  €  Hp(12),  i  =  1,2  and 


(4.12a) 


II v  II  <  Ctk_1llult 

1  Hx(n)  Hk(n) 


(4.12b) 


2  k-2 

II  v  H  „  <  Ct  Hull 

Hz(n)  Hk(n) 


with  C  independent  of  u.  Then  by  Theorem  4.1  there  exists  zp  €  Sp  such 
that 


Hz  -v„fl  ,  <  C p  ^  II v 

P  2  H1  (12) 


-1  k-2 

A  0  <  Cp  t  Hull 

H  (12)  H  (12) 
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-V-i-.'VVV' 


V  nV  A.  /«  aV  A  V-  l^kA  M*-  . 


Choosing  t  *  1/p,  (4.12a)  gives 


8  v  0 
1  1 
H1  (ft) 


Cp-(k-l  )lijg 


Hk(ft) 


We  get  by  the  triangle  inequality 


(4.13)  It z  -ut  <  Cp”^k“1^tull 

P  H*(ft)  Hk(ft) 

i.e.  the  estimate  (4.1). 

The  assumed  splitting  (4.11),  (4.12)  is  equivalent  to  the  definition  of 

1  O  1 

the  interpolated  space  (HD(ft),H  (ft)  fl  H^(ft))  defined  by  the  K-raethod 

9  ,® 

(see  [6]). 

We  have  shown  in  [3]  that 


(H*(ft),H2(ft)  n  Hp(ft)) 

9  ,  °° 

*  (H1(ft),H2(ft))0  w  n  ^(ft) 


-  Bzl(n)  n  Hi(n)  - 


-  (H  1  (ft)  ,H  2(ft)  ) 0  2n  Hp(ft) 


-  Hk(ft) R  (ft) 

1+0  k 

where  B2  ^  is  the  usual  Besov  space.  Hence  (4.11),  (4.12)  hold  if  u  £  H^( 

and  Theorem  4.2  is  proven. 

Remark  4. 1 .  We  mention  that  we  have  proven  slightly  more  then  the  claim  of 

Theorem  4.1,  namely  that  we  can  use  the  Besov  space  B_  ^(ft)  instead  of 

*•  * 

H  (ft) .  This  of  course  follows  easily  from  the  interpolation  theorv. 


Remark  4.2.  The  assumption  ok  <  2ir  (i.e.  the  exclusion  of  the  slit  domain) 
was  made  here  only  because  it  is  used  in  a  result  from  [3]  quoted  by  us. 

4.3.  The  rate  of  convergence  of  the  p-version  of  the  finite  element  method 

THEOREM  4.3.  Let  u  €  ,  k  >  1  be  the  solution  of  (2.2),  (2.3).  Assume 

further  that  g  is  such  that 

u  =  u1  +  u2 
k,  k2 

uL  €  H  £<Q),  u2€H  z(n),  k2  >  3/2 

and  that  ft  is  a  Lipschitz  domain  if  k^  <  3/2.  Let  up  be  the  finite  element 
solution  based  on  the  p-version  defined  in  Section  2.  Then 

(4.14)  lu-u  II  ,  <  Cp_(k~1)llull  ,  k  *  min(k,,k2). 

P  Hx(n)  H  (a) 

The  constant  C  is  independent  of  p  and  u.  It  depends  on  the  factorization 
of  fi  into  elements  0^. 

Proof .  If  g  ■  0,  then  (4.14)  follows  immediately  from  Theorem  4.2  because 

llu_u  II  .  <  Cllu-z  H  ,  . 

P  H  1  (12)  P  H^Jl) 

If  g  ^  0  then  denote  by  Up  the  exact  solution  of  the  problem  (2.2),  (2.3) 
when  replacing  g  by  gp.  Denoting  w  =  u  -  Up ,  function  w  obviously 
satisfies 


-Aw  +  w  =  0 


3w 

3n 


0  on  rN 
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.V.*h  -V.V  .V  a  -%'Ov  .V.v'aIii.V^  .n  a"a  .Va’a'a’a'a'UWa  -v*.v  a1 


A  .v 


,*v-s 


(I) 


V 


x 

r 


I 


B 


g  -  gp  on  r  . 


Using  Lemma  3.2  and  the  same  argument  as  in  Theorem  4.1,  i.e.  the  extension  by 


the  function  $p  in  (4.4)  we  conclude  that  there  is  a  function  w  €  H*(ft)  such 


that 


-(k-1  ) 

Uwa  <  Cp  v  '»u« 

ir(n)  Hk(ft) 


and 


w  *  u  «  g  -  gp  on  T  . 


Hence 


(4.15) 


-(k-1 ) 

Hull  <  Cp  v  'Hull 

Hx(n)  Hk(n) 


because  w  minimizes  H*H  .  among  all  functions  with  trace  g  -  gn  on 

n  (ft)  P 


Tu.  By  Theorem  4.2  and  a  basic  property  of  the  Finite  Element  Method,  we  have 


Hu  -U  H  ,  <  C II z  -U  II  . 

p  PrV)  P  p  u  1 


H  (ft) 


<  C  ( II  z  ~u  II  .  +  H  u-U  II  . 

P  H  (ft)  P  H  (ft) 


<  Cp 


-(k-1) 


H  (ft) 


and  Theorem  4.3  is  proven. 
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5.  THE  CONVERGENCE  RATE  OF  THE  p- VERS ION:  THE  CASE  OF  THE  SINGULAR  SOLUTION 

In  Section  4  we  analyzed  the  rate  of  convergence  of  the  p-version  under 
the  assumption  that  u^  *  0  in  (2.4).  In  this  section  we  will  consider  in 
detail  the  case  when  u  =  u3^* 


5.1.  An  approximation  result 

Let  Q  =  (-1,1)  x  (-1,1)  as  in  Section  3.  Denote  3L  =  x,  +  1 ,  i  * 
1,2  and  let  for  k>1,  0  <  p  <  1, 


S 

< 


S 


P 

K 


% 


R 


K 


< 


{x  €  Q  I  i  Sj  <x2  <  kXj} 

s^  n  {x  |  +  x2  ^  p2} 

{x  I  0  <  *  <  1,  0  <  3?2  <  1} 

{x  I  0  <  X1  <  V2  ,  o  <  5?2  <  l/2  } 

S<  n  % 

SK  nq0 

sK  n  Cf0* 


Let  <n  >  <  >  1.  Fig.  5.1  shows  the  domains  under  consideration. 


Fig.  5.1 


2  ~2  ~2 

Let  (r,0)  be  the  polar  coordinates  with  origin  (-1,-1);  r  =  +  *2 ’ 

x2 

0  =  arctan( — )  . 

*1 

Let 

(5.1)  5(jKj  ,X2)  -  (Xl-icX2)(idt1-S£2)  -  r2»jO). 


Obviously  <t>  ^  ( 0 )  is  an  analytic  function  in  9,  and  £  is  a  polynomial  which 
vanishes  on  the  lines  Ki  =  <X£  and  Xj  =  —  ^2* 


Let,  for  a  >  0,  y  >  0 


(5.2) 


uCx^.S^)  *  ra|log  r|Y  x(r)i>(9) 


where  #(9),  x(r)  are  sufficiently  smooth  functions  (e.g.  C*  functions)  and 


x(r)  -  1  for  0  <  r  <  4 


X(r) 


,  2p  n  1 

for  -j-  <  r,  0  <  p  <  — 


is  a  function  defined  on  Q.  We  shall  assume  that  u  vanishes  on  the  lines 
3t,  M  xK-  and  H ,  «  —  and  has  support  in  .  Then 

1  l  1  <4 


(5.3) 


Uq ( x  ^ jx 2  ^ 


u(2^ , Sf2 ) 
5 (x^  «x2 ^ 


ra  "|  log  r|Yx(r)'^(0) 


where  ^(G)  is  once  more  smooth  (e.g.  C*  function). 
Now  we  can  write 


(5.4) 


u(x1,x2)  -  C(x1,x2)u0(x1,x2). 


The  main  result  of  this  section  is 


THEOREM  5.1.  Let  u  be  given  by  (5.2).  Then  there  exists  zn  £  P p+->(Q)  such 


1 


that  z„  *  0  on  the  lines  3 1,  *  and  xL  *  —  3L  ,  and  for  <n  >  <, 

P  12  1  <2  0 


(5.5) 


lu-z  ll 

p  Hla  ) 

*0 


<  c|  log  p|  Yp 


where  C  is  a  constant  independent  of  p. 

We  will  need  a  series  of  lemmas  to  prove  Theorem  5.1 
Let  u>(r),  0  <  r  <  «®  be  a  C*  function  satisfving 


30 


<i)(r) 


0 


for  0  <  r  <  1 


co ( r )  *  1  for  2  <  r  <  ®. 


Further,  for  any  A  >  0  let 


(5.6) 

and 

(5.7a) 


u>A(r) 


»(-£) 

A 


A 

0)  ur 


(5.7b) 


w  *  (1-0)  )u. 


Then  obviously 
(5.8) 


u0 


v  +  w. 


It  can  be  readily  seen  that 

v  =  0  for  0  <  r  <  A 

w  -  0  for  r  >  2A. 


LEMMA  5.1.  Let  k  =  k  j 
exists  a  constant  C(k) 

ak 

i  3  v 
1  k2 

(5.9)  3xl  3x2 


+  k£  where  0  <  kj,  k2  <  k 
such  that  for  x  =  (xj.xj)  ? 

<  C(k)|log  A| Y( l+xi)a  “  ^ 

»  0 


are  integers.  Then  there 


Proof .  We  have 

1 1— -|  <  C[|r3  3|log  r  |  Yx(  r)4>(0)wA(  r)  i  +  |r3  ^|log  r|Y  *  x(  r)  *(  0)o)A(  r ) 


I  m 


+  j  ra_2  |  log  r|Yx'(r)i|»(9)uA(r)|  +  |ra'2|log  r  |  Yx(  r)«K©)</' (  r )  -|  ]  . 


Note  that  the  third  term  on  the  right  hand  side  is  0  except  for 
<  r  <  and  the  fourth  term  is  zero  for  r  >  2A.  Hence  for  A  <  y  we  have 

I 111  _  o  for  0  <  r  <  A 

'ar1 

<  C | log  A|Yra-3  for  A  <  r  <  2A 

<  C|log  a|Y  max(ra  2,r°  3)  for  r  >  2A,  y  <  <  ■ 

<  C|log  A|Yra  3  for  r  >  2A,  r  ^  (-y,  -j— / 


Hence 

||^|  <  C | log  A | Y r a  3 

for  all  r. 

This  process  can  be  repeated  to  obtain 


(5.10) 


|-~|  <  C(k)  |  log  A|Yr°'2'k. 

3r 


In  S  we  have 


D^(l+Xi)  <  r  <■  (1+x^)D2 


3  r 

3xl 

cos 

0, 

3  r 
3x2 

QJ 

CD 

1 

sin 

_e 
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3x  ^ 

r 

> 

rJ 

X 

ro 

sin  0 


cos  0 


and  ip (G )  is  smooth.  Hence  (5.10)  gives  immediately  (5.9) 


In  what  follows  we  will  assume  that  v  satisfies  (5.9)  and  not  the 
explicit  form  (5.3),  (5.7a). 

Let 


(5.11) 


v(Xj  ,x2) 


l  l  a  P  (x  )P  (x2) 
i=0  j-0  1  J 


where  P^(x£)  =*  Pi(x£,B,6),  6  >  -  y  are  Jacobi  polynomials  of  index  B  which 

will  be  determined  later.  Then 

aij  =  CiCj(i+l)(j  +  l)  /  /  v(x1,x2)Pi(x1)P.(x2)(l-xf)6(l-x2)8dx1dx2 

where  are  bounded  from  above  and  below  independently  of  i,  j  but 

depending  on  6  (see  [12],  p.  841,  formula  7.391.1). 

Define 


(5.12) 


P  P 

vp(*l»x2)  a  aijPt(x1)Pj(x2) 


(5.13a) 


(x2)  .  I  .ijPjfxj), 


j-0 


(5.13b) 


with 


(5.14) 


b[pl(x„) 
l  2 


l  aijPj(x2) 

j-0 


+  1 


bi(x2)  =  Ci(i+1)  /  v(xj ,x2)( l-x^)SPi(x1 )dx1 


It  can  be  readily  seen  that 


(5.15a) 


v  -  I  b1(x2)Pi(x1) 
i-0 
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£ 

5$ 

i 


— -  has  support  in  and  is  zero  for  Xj  lying  outside  the  interval 


Ij(x2)  =  [-1  +  —  (l+x2),  -1  +  <q  (l+x2)]. 


Using  Lemma  5.1  we  get 


* 


i 


d^  (x~) 
1  1 


C(i+1)^2  /  |log  a|Y(1+x  ) 


B  _  I 

Cl  C.  111  *  A  / 

2  4 


1  ^(x2) 


r\  .8  1  ,  . 

l/„,  v  a-2-m+j-^+l 

<  Ci'2  I  log  A| Y(l+x  ) 


8  5 


provided  that  ct  —  m  +  —  —  <  0  which  is  (5.19)  for  i  >  1.  The  case  i  -  0 

is  verified  separately  (using  the  fact  that  8  >  -  V2  ) • 


LEMMA  5.4.  Let  a  -  m  +  j 


-7  <  0.  Then 
4 


(5.20) 


d™!) .  (  x  «  ) 
1  1  i 


„  .tY  a-m+  -r 

-g  -i-  (1+xJ  2 

th  2 


(i+1) 


Proof .  We  have  (see  [12],  p.  1039,  Formula  8.964) 


(5.21) 


(l-xpP'(x  )  -  2x1(B+l)P'(x1)  +  i(i+28+l  )P.(x1)  =  0, 


Multiplying  (5.21)  by  (l-x[)B  we  get 


(5.22) 


-i(i+2B+l)(l-x2)3P  (x) 


(ci-x2)1+3p:(x  )) 
dx  1  1  1 


Hence,  differentiating  (5.14)  m  times,  using  (5.22)  and  integrating  by  parts  we 


(5.23) 


dmb. (x_) 
1  1  2 


+  1  „m 


‘  c  (wry  U,  ^7^7  ((|-x?)  plui))dxi 


-N  h  -  ■  -  *  *  *  -  ■  »  p  -  *  *•  •*  k.V  -VJV -V  »  '  ^  .V  *  r  .VA.V.  s,~»VvV -V 


ijt/  l*r(r-T=(l-N>  )' 


-1  “1  3x  3x 


VI 


As  in  the  proof  of  the  previous  lemma,  the  integrand  is  2ero  outside  of  i^x2 
Further, 


(5.24) 


,  ,  m+1  . 

(- - —  (l-xf) 

9xi 


1+8 


_  1+0  ,m+2  nrt-1 

<  *  |2(i+e>  s  v 


9x^  9x^ 


3x^  3x^ 


B  ^'V  V 


<  C(m)( l+x1)a-m+®  3| log  A|Y 


where  we  used  Lemma  5.1  and  the  fact  that  (l-x^)  <  2(l+x^)  on  I^(x2) 
by  Lemma  5.2,  (5.23),  (5.24)  we  get 


(5.25) 


dmbi ( x2  ) 

dx” 


<  c  Ll-°-g-34i  j'  (1+X1} 

(i+l)J//  IL(x2) 


^  8  13 
o-m+  -=»  -  — 

2  “4 


dx . 


Hence 


£li°ij^(i+x 

(i+D3/2  (  2 


3  9 


provided  that  a-m  +  Y~-^-<0.  Combining  (5.19)  and  (5. 25) , 


(5.26) 


.  B  5 

,  \  n ®  nn-  o  a  1  / 

1  1  <  C|log  a|Y(1+x7)  min{(l+l) 


d^ .  (x„) 


(1+x  J 


dx. 


( i+1 ) 


3/2' 


Using  the  logarithmic  inequality 


min{a  ,b}  <  a ^  b^ 


yields  (5.20). 


given  in  (5.17).  We  have 


Let  us  analyze  now 


v 

P 


Because 


Pp(xi>x2) 


^  [ bi(x2)_b|P^ (x2> ]P ^(x ^ ) 


f  [bi(x2)~b|p^ (x2) ]P'(x1 ) . 


P'(x,0 ,0) 


j  ( 20+i+l )Pi_1 (x , 0+1 ,0+1 ) 


(see  [12],  p.  895  formula  8.961.4),  we  obtain  for  0  <  m  <  p  +  1 


(5.27) 


+1  +1  3p  2  ?  6+1  9  3 

^  =  J  (/  (xl  »X2^  (l_xl'j  dxi)(1_x2)  dx2 


+1  D  r  I  2  p  S 

<  C  /  (  l  i(bi(x2)-b^p  (x2))  (l-x2)  dx2 


-1  i-1 


<  C  f  i  /  (  l  a  P  (x  ))  (1-xJ)  dx 
l-l  -1  p+1  133 


00  a,  . 


c  1  I  1  -f  ^  h  I 

i=l  j=p+l 


CO  a  ( j+m) ! 


-LI 


j  p2®  it,  *  j-p-n 


p  +1  dmb,(x7) 

1  v  r  r  1  - 


0+m 


<  c  -s  l  1  /  ( — — )  (l"xP  dx; 

p  i=l  -1  dx™ 


Using  (5.14)  we  see  that  the  support  of  b i ( x 2)  lies  in  It  =  [-1+4  sin 

1  -  3  5/ 

where  tan  0n  =  — .  Hence  from  (5.27)  and  (5.20)  for  a  -  m  +  —  -  —  <0 
n  <  - 

0 
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a 


p  i-1  -l+A  sin  0q  dX2 


log  A 


i=l  -l+A  sin  0 


U  , 

/  t  n+x  ) 


2(a-m+ 


I .  .  |  2y  2a-m+20-  ■=■ 

<  C  I10?  fl  i  2 

2m- 1 
P 


provided  that  2a  -  m  +  23  — ^  <  0.  Analogously,  using  (5.19)  for 

o  5 

a-m  +  y-  -£-<0  and  0  <  m  <  p  +  1 , 


+1+1  ~  -  8  «  6 
(5.29)  A2  =  /  (/  (Pp)  (1-xJ  dx1)(l-x2)  dx2 


+1  P  ,  f  i  2  _  8 

<  C  /  I  —  (bi(x2)  -  b| P  ( *2 ) )  (1_x2)  ^ 


r  P  ,  0  d^b^x,)  2  _  8+m 

-S  i  T  /  C — t1”1  >  d* 

pi-1  -l+A  sin  0Q  dx“ 


1  .../ 


„  2(a-m+  4-4) 

i|log  A|  T(1+xJ  ( 1+x , 


p  i-1  -l+A  sin  0 


„  2a-ra+2B- 

<  “2^T  |log  ^ 

P 


provided  that  2a  -  m  +  23  -  <  0. 

Similarly  we  estimate  the  term  a  in  (5.17)  with  0  <  m  <  p  +  1 


0p(xl,x2) 


.  *  bi(x2)Pi(x2) 
i=p+l 


(5.30) 


+1  +1  3a  ,  8+1  ,  8 

-l  ^1  ^  ^2 


+1  ®  .  »  8 

r  r  r  ^/_.  \\r ,  ^ 


C  /  (  l  ib  (x  ))(l-x  )  dx 

-1  i-p+1 


r  +1  »  b. (x_ ) ( i+ra) !  „  6  . 

<  “XT“  TT"  f  y  -1,.2  ... -  fl-xll  dS2 


2(m-n  L,  <i-«m  (1_x2' 

p  -1  i*p+l 


_  +1  +1  -m  2  .  S-Tn  .8 

<  ~2'(  m-T )'  /  U  Hi)  (1_xl)  ^l)^"^)  dx2* 

P  -1  “I  OXj 


Since  the  support  of  v  lies  in  R,.  -  S  we  can  use  Lemma  5.1  and  obtain 

*0  *0 


with  1^  *  [-1  +  —  (l+x2),  -1  +  <q(1+x2)] 


(5.31) 


,  „  | log  a|2y  ,°  ,  ,,.  ,2<«-2-«),.  2B+m  ,,  2?  .  , 

B1  ‘  C  2(S-1  J,  .  {  (l-x2)  dx  dx 

p  -1+A  sin0Q  I1 


rf1  /°  u+x,)2“-2s-3dx: 


2(m-D  ...  a  Vi  2 

p  -1+A  sin  0n 


„  „  I  log  Al  Y  . 2a-m+28"2 

*  c  2(m-rr" A 

p 


provided  that  2a  -  m  +  2p  -  2  <  0. 


Similarly  for  0  <  ra  <  p  +  l 


+  1  +1 


(5.32) 


2,-.  2s3,.  2x 3 


/  /  Cl"xi)  ^1_ 


x- 1  dx, dx 


-1  -l 


2  i  ^  i ? 


vm wwm: 


_  +1  +1  »m  2  0  8+m  0  8 

<  s!  I,  (sf)  *‘ld*2 

p  “1  “1  1 


<  C  |log  s|2?  42a-»+2B-2 


provided  that  2a-m+28~2<0. 

We  will  summarize  (5.29)  -  (5.32). 


LEMMA  5.5.  Let  pp  and  ap  be  defined  in  (5.17).  Then  for  0  <  m  <  p  +  1 
and  k  *  0 , 1 


(5.33) 


+1  +1  3  p  2  0  8-Hc  ,  8 

/  /  (— J2)  (1_x2)  dx1dx2 

-1  -1  3x* 


i .  ,  |  2y  2a-m+28-  -  k 

<  cJJfiL*  2 

2nr-l 


8  5  3 

provided  that  a-m  +  y-  -£-<0  and  2a  -  m  +  28  _  y  _  k  <  0 


(5.34) 


+1  +1  3a  2  _  6+k  ,8 

/  /  (— JJ2)  (1_x2)  dx1dx2 

-1  -1  3x^ 


| log  A [ 2y  2a-m+2B-2 

L  2(m-k)  A 
P 


provided  that  2a  -  m  +  28  -  2  <  0.  The  constant  C  is  independent  of  A, 


p  but  depends  on  a,  8,  Yt  ra> 


ra  -  r  n  q^a 

K  K  0 


where 


*2  >  “1  +  2 A } 


,2A 


{x^  >  -1  +  2A, 


Further,  for  f(xj,x2),  xlfx2  €  Q 

(5.35) 

fA^Xl’X2^  *  f^X 

and  f^(xlt 

x2)  -  0  on  Q  -  QqA. 

LEMMA  5.6. 

Let  £(x^,x2)  be  given 

(5.36) 

1  »x2) 

HA 

(5.37) 

Proof.  For 

(xl5x2)  €  R 

f  (x^**2A  ,x2~2A)  ,  (xltx2)  6  Q" 


2a 


1  x2+1 

K0  <  xl+1  <  <0 


and 


for  (x1,x2)  ^  R^ 


2A  <  1  +  Xj . 


Hence 

(4.38) 


Ua(x1x2)I  * 

“  | { 1+x1-2A-k(1+x2-2A))(<(1+x1-2A)  -  (l+x2~2A))| 

3  Id+Xj)  -  <(l+x2)  +  2A(<-1  )  I  I  <(l+x1 )  -  (l+x2)  +  2 A ( 1  —  <  )  | 
<  (|<Q(l+x2)|  +  |k(1+x2)|  +  |  ( tc-1  )<0(i+x2)  |  ) 
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X  (|k(1+Xi)|  +  iKgd+x^l  +  ICie-Dd+Xj)  I  J 


<  (1+X2>(  1+x^  )(<Q+tc+(<-l  )<q)(<+k^+ic-1  ) . 


Because  on  R„ 


(l+x2)  <  C(l-x2> 


(1+x^  <i  C(l-xp 


(5.38)  yields  immediately  (5.36).  (5.37)  can  be  proven  in  an  analogous  way, 


LEMMA  5.7.  Let  v  satisfy  (5.9)  and  vp  be  given  by  (5.12).  Then  for  A 


(5.39) 


(v-v  )ll  . 

4  p  hV  ) 

0 


<  C j log  p|Yp  2a 


where  C  is  independent  of  p. 


Proof.  As  in  (5.17)  v  -  v  *  a  +  p  .  Let  us  estimate  first 

II5.P  II  ,  .  .  To  this  end  we  estimate  H5.P  H  .  .  It —  P  i!  . 

4  P  HV  )  4  P  L,(R4  )  9X1  P  L,(R4  ) 

*0  2  *0  2  *0 


P  3 

R  5 .  t— HI  .  The  estimates  of  the  terras  involving  r —  follow  from  the 

A  dx,  r  A  .  0  9x„ 


1  L2(K  }* 
0 


symmetry  of  Xj  and  X2, 


Using  (5.36)  we  get 


(5.40) 


3p  2 


3p  2 


^A  3xi  l2(ra  ) 
*0 


5  A  ^  dx  1  dx 2 
R A  1 

Ko 


2  2  3  p  2 

<  c  //  (l"xi)  (l-x2)  dXjdx2. 

RA  1 

K0 


f»  ^  **.  V  •.,f  •'  *  *•'•''*  *  „  , 


■'CiuVj.'.  tS.\S  kS 


Because  on  R  , 

^  r\ 


0  <  C  < 


2A  ’ 


using  Lemma  5.6  we  get  for  S  >  2 


,  s+i  2  8 

+1  +1  ( 1— x  )  (l-x_)  3p  2 

»!«<=/  /,  -d-  Csf)  *^2 

-1  -1  A  A  I 


0  2a-m+2B-  4  -26+3 

< 

p 


C  I  log  A  I 2y  - 


2a-(m-  V2  ) 
2(m-  V2  ) 


8  5  5 

provided  that  a-m  +  y-  ^-<0  and  2a-m+2B-y<0. 

-2 

Choosing  m  large  enough  and  A  *  p  we  get 


(5.41) 


Similarly 


Dr  <  C| log  p| Zyp  4a. 


3Ca  2 

D?  “  p  J  A 

2  3X1  P  L,(R&  ) 

2  <0 


<  C  /  /  (1-x2)  (pp)2  dx1dx2 
R 

*0 


8  0 

+1  *1  (l-x?)  (l-Xj)  2 

‘  C  V,  - <°P>  dXldX2 

-1  -1  A  A 


f  mV  mV  -h\  Wj*S  »  *-M  *_■***  7^VfcV_»  R-  *V  »N  AjV  *V  A  i.V^ 


C  I  log  A 


8  5 


2y  A 


^2a-m+  V2 
2  (nr-  V2  ) 


provided  that  a-m  +  y-  -^-<0  and  2a  -  m  +  28 


7  <  °- 


Hence  for  m 


large  enough 


(5.42) 


(5.43) 


D2  <  C |  log  p|2V4a- 


B  B 

+1  +1  (1-x2)  (1-x2) 

*  <  /  /  - 2fi -  PndX]dX2 

(ra  )  -1  -1  A2e  P  1  2 

K0 


<  C|log  p|2V4a'4. 


Combining  (5.41),  (5.42),  (5.43)  we  get 


(5.44) 


"5aV  1  A  <  c  I  lo8  P|YP  2  • 

A  P  hx(ra  ) 

K0 


The  estimate  for  II £,a  H  .  .  can  be  obtained  quite  analogouslv. 

1  P  bV  ) 

K0 

Let  now  a  >  1  not  be  an  integer  and  k  =  [a]  be  the  largest  integer 


less  than  a.  For  0  <  x  <  k,  r  integer,  let  v^T- 


denote  the 


derivative  of  v  (see  5.7a)  along  the  direction  n  where  n  is  the  unit 
vector  along  the  line  x2  ”  Xj .  Then  we  see  that  v^T^  satisfies  Lemma  5.1 
with  a  replaced  by  a  -  t  >  0.  Hence  using  Lemma  5.7  we  get 


(5.45) 


r  -  »r  w 


r  .  t.t  «-  y~.  S'  v-  J-  «r  «r  j-  ■ 


■S4(v[Tl-v'Tl)i  ,  .  <  C|  log  p|V2(“-T)- 


P  H 1  (R  ) 
*0 


Let  tu“  be  defined  by  (5.6)  and  be  its  translation  given  by  (5.35) 


Then  (see  5.8) 


(5.46) 


U0A  *  U0AUlA  +  U0A(1_UiA) 


v  +  v  . 

A  A 


Because  u  €  H^. ( R  ),  then  u.  €  h!(R  )  and  hence 
0  *0  A  0  *0 


CAWA 


U.(l-ui^)  f  H  1  (RA  ) 

A  A  <Q 


.-2 


LEMMA  5.8.  Let  A  =  P  ,  2  *  2/2  A.  Then  for  k  =  [a],  a  >  0  noninteger 


(5.47) 


I  (-D^v'1'). 

i=0  H 1  (RA  ) 

*0 


<  C | log  p|Tp  2“ 


(5.48) 


«5awa«  ,  A  <  c|log  p|Yp 

A  hx(ra  ) 

0 


where  C  does  not  depend  on  p,  A, 


Proof .  By  Taylor's  theorem,  for  any  (x^,X2)  €  R  and  s  =  0,1,  we  get 


using  Lemma  5.1 


k  ,1  41  .ID 


-sK  (-olfr 


,x. 


3x 


i-0 


,  , ,  .  .  k+1  -s 


45 


<  CA^Cl+x^-^llog  a|Y 


where  0  <  1 0 j  <  2A.  Hence  using  Lemma  5.6  we  get  for  A  =  p 


i  Ji]^2 


w'Vhlvj. 


L"(R°  ) 
<0 


<  C  //  (1-Xj)  (l-X2)2SA2(k+1)|log  Al^Cl+x^ 

R* 

Ko 


2ct-6-2s-2k 


<  CA 


2(k+1>|log4|2’'  /  (l+x,)2(“-1‘-2)+1dx, 

2A 


<  C|log  A|2YA2(k+1)+2^a“k-1) 

<  c |  log  pi  2yp  4a. 


In  the  above  inequality  we  used  the  obvious  fact  that  a  -  k  -  1  <0.  The 
other  terms  in  (5.47)  can  be  bounded  analogously  and  the  first  part  of  Lemma 
5.8  is  proven. 

Let  us  now  prove  (5.48).  It  is  easy  to  see  that 


"W  1,  A  , 
H  (R  ) 

<0 


H  (T) 


where  T  =  {(r,0)j,  0  <  r  <  2A,  0  <0  <  -j}  .  Using  (5.1)  we  have 


(5.49) 


I  £(  r  ,0)  |  <  Cr^ 


(5.50) 


iff-  (r, 0)|  <  Cr- 


Further,  by  (5.3) 


completely  proven. 


Now  we  can  prove  our  main  result 


Proof  of  Theorem  5.1 


Let  S  be  the  translation  of  obtained  by  the  transformation 


*  xi  -  2A.  Let 


i!  p  J’ 


k  =  [a] . 


Then  zpA  €  Pp+2(Q)  and  zpA  “  q  on  c^e  sides  of  S<.  We  have 

k  .  A  .  r  ,  I 

Hu  -z  .11  ,  .  =  »£A(unA  -  T  (-1)  tt  vl  J)n  .  , 

A  P^H1(RA)  ^ A^  OA  ^  i!  p  ^  H 1  (RA) 

*0  Ko 


'C.w.H  ,  .  +  US  A(v 

1  4  H  V  )  4  4 

‘0 


k  .  A .  r  . , 

l  (-n1^  »Illh  ,  4 

i=0  H 1  (R  ) 

Ko 


♦  I  !  4 

i»0  p  H  (R  ) 


by  (5.45),  (5.47)  and  (5.48) 


<  C | log  p|Tp  2a. 


Translating  back  to  SK  we  get  the  theorem. 

Remark  5.1.  We  have  proven  more  than  Theorem  5.1  claims.  It  is  sufficient  to 
assume  that  v  and  w  defined  by  (5.7a),  (5.7b)  satisfy  (5.9)  and  (5.51), 
(5.52)  respectively. 

Remark  5,2.  It  is  easy  to  see  from  the  proof  that  the  internal  angle 

of  y ^  and  could  be  equal  to  2n ,  i.e.  that  we  may  also  consider  the 

slit  domain. 


[i] 


5.2  The  rate  of  convergence 

We  return  now  to  the  problem  of  approximation  of  the  function  u^ 

(2.4). 

Let  the  vertex  be  at  the  origin  0.  The  part  of  the  domain  £ 

containing  the  elements  with  vertices  at  0  is  shown  in  Fig.  5.2. 


We  will  assume  that  the  lines  OB.  have  the  coordinate  0^  and  that 

3  J 

-  D  -  N 

OBj  c  r  ,  c  r  .  Although  we  assume  that  we  have  only  triangular 

elements,  the  case  when  elements  are  parallelograms  does  not  change  the 
argument . 


m 

ra 

Let  £2 

=  U 
i=l 

V 

i=2  1  1+1 

*  r . 

Denote 

Op  -  <*  1 

and  assume  that  V 

p0 

c  £2, 

0  <  pQ  <  l. 

Now 

we 

prove 

THEOREM  5.2. 

Let  u 

be 

the  function 

given 

by 

(5.2) 

with  p  < 

suf  f iciently 

small . 

Then 

there  exists 

Zp  € 

H1 

(£2) 

satisfying 

1  *  l ,  . . . ,m, 

ZP  =  ° 

on 

OB}  and  T 

and 

'P  P 


(5.53) 


I  u-z  II  <  C  |  log  p  |  ^p  “a 

P  H  1  (£2 ) 


where  C  is  independent  of  p. 


1.  Assume  first  that  $(0^)  “0,  j  =*  l,...,m+l.  Denote  $^(0)  to  b 
an  extension  of  $(©)  onto  (®j“A,0^  +  ^+A)  where  |2A+0^.+^-0j!  <  tt  and 
fj(0)  =  0  on  (0j~A »© j”  f)  and  (0^+  j,0..+A) 


0j-A 


0  V  SP 


eM+A 


Fig.  5.3 


Let  S  =  { (r ,0) | Oj  —A  <  0  <  0j+1+A}  and  S  =  { ( r , 0 ) | 3 ^  <  0  <  0j+]}. 

Denote  by  Uj  the  function  given  by  (5.2)  when  $(©)  is  replaced  by  $^(3) 

and  extend  uj  by  zero.  Let  now  T  be  the  linear  mapping  which  maps  S 

onto  R  and  S  onto  R  .  Denote  T.  =  T(T-)  and  assume  that 

*0  <  3  J 

T.  C  R  CQ 
3  <  0 

~  -s* 

The  mapping  T  transforms  u4  into  u.  on  T..  Denoting  bv  n.  the 

J  J  3  ~  3 

_ _  u . 

linear  function  which  is  zero  on  T(b.B..,),  the  function  — J-  satisfies 

3  3+l  n. 

obviously  the  conditions  mentioned  in  the  Remark  5.1  to  Theorem  5.1. 


Therefore , 


can  be  approximated  by  a  function  z  satisfying  (5.53)  on 


satisfies  (5.53)  too.  Hence  (5.53)  is  proven  in 


^  -k 

T.  and  hence  z  h  .  =  z  ,  , 
j  P  J  P+1 

the  case  that  $(0j)  =0,  j  =  l,...,m. 


2.  Now  we  will  consider  the  case  when  $(6^)  =  0  for  .1  t  j q  . 
the  elements  T,  as  shown  in  Fig.  5.4 

Jo  1  Jo 


Consider 


If  the  angle  |0.  . -0 .  .  |  <  tt  ,  then  we  can  proceed  exactly  in  the 

J0  J0 

same  way  as  before  onlv  replacing  n.  bv  n.  ,n.  .  Hence  we  have  to 

J  '  J0  30 

consider  the  case  when  0 


.  .  -  0  >  *  (see  Fig.  5.5) 

J0  J0-l 


6.  THE  CONVERGENCE  OF  THE  p-VERSION  OF  THE  FINITE  ELEMENT  METHOD 


In  this  section  we  will  summarize  the  results  we  have  proven  and  further 
generalize  them. 

6.1.  The  case  with  triangular  and  parallelogram  elements 
We  have 


THEOREM  6.1.  Let  ft  be  the  polygonal  domain  as  introduced  in  Section  2. 
Consider  the  problem  (2.2)  (2.3)  and  assume  that  the  solution  can  be  written 
in  the  form  (2.4)  (2.5)  with  k  >  3/2  and  that  ft  is  Lipschitz  for  q  <  3/2. 


Assume  that  Up  is  the  finite  element  solution  as  described  in  Section 


2.3  with  triangular  or  parallelogram  elements.  Then 


(6.1) 


iu-u 


”  h‘w) 


<  Cp  U | log  p| 


where 


(6.2a) 


=  min(q-l,  k-1,  2a[*^)  =  min(q -1  ,k-l ,  2a  ) 

i  1  1 


v  = 


max  Y  if  U 


=  2a 


lj] 


(6.2b) 


0  otherwise 


(6.2c) 


lui:l  q  +  "u2"  k  +  I  lC£ 

Hq(ft)  H  (ft)  i  ,£ 


l  i) 


The  theorem  follows  immediately  from  Theorem  4.3  and  5.2. 

Although  we  have^dealt  with  the  model  problem  (2.2),  (2.3)  only,  it  is 
obvious  that  the  theorem  holds  for  any  second  order  elliptic  problem  if  the 


53 


solution  has  the  fora  (2.4),  (2.5)  or  when  (2.5)  is  different  but  has  the  same 
character  concerning  the  growth  of  its  derivatives. 

6.2.  The  case  of  curved  elements 

So  far  we  assumed  that  the  elements  are  triangular  or  parallelogram. 

The  obtained  results  can  be  immediately  generalized  to  the  case  of  curvilinear 
triangles  and  quadrilaterals  which  can  be  mapped  individually  on  the  standard 
triangle  or  square  by  a  mapping  which  is  one-to-one  and  sufficiently  smooth 
(in  both  directions).  (In  practice  this  is  always  achieved.) 

We  proceed  then  in  the  usual  manner  by  approximating  the  function  on 
standard  squares  and  triangles. 

Theorem  6.1  holds  then  without  changes.  The  function  u^  is  of  course 
defined  now  on  the  straight  line  triangles  (squares)  and  does  not  have  the 
explicit  form  given  in  (2.5)  but  possesses  the  same  character. 

So  far  we  assumed  that  all  elements  are  of  the  same  degree.  The  modifi¬ 
cations  in  our  proofs  and  results  when  the  degrees  are  different  in  different 
elements  are  obvious. 
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